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Abstract 
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1 Introduction 



For integrable models, the relationship between classical and quantum systems is in no 
way exhausted by their correspondence in the classical limit. As is now well known, clas- 
sical integrable equations often appear in quantum integrable problems as exact relations 
even for h ^ 0. 

One important example of this general phenomenon was investigated in [U|2l[3], where 
it was shown that the spectrum of commuting transfer matrices (integrals of motion) in 
quantum integrable models can be found in terms of discrete classical dynamics, also 
integrable, defined in the space whose points label the commuting transfer matrices. 
For integrable GL(if )-invariant spin chains coordinates in this space are parameters 
specifying finite-dimensional irreducible representations of the group GL(K) and the 
spectral parameter. The classical dynamics in this space is generated by functional 
relations for the transfer matrices established in [U [5j E] for the ordinary bosonic case 
and extended to the supersymmetric case in [Jj. Among them the most important is 
the bilinear functional equation for the eigenvalues of the transfer matrices (T-functions) 
which has the form of the Hirota bilinear difference equation [8]. For brevity, we call it 
the TT-relation. It is the starting point of our approach. 

The Hirota equation is probably the most famous equation in the theory of classi- 
cal integrable systems on the lattice. It provides a universal integrable discretization of 
various soliton equations and, at the same time, is a generating equation for their hierar- 
chies. It is involved in a large body of integrable problems, classical and quantum. Like 
other soliton equations, the Hirota equation admits (auto) Backlund transformations, 
i.e., transformations that send any solution to another solution of the same equation. 
They allow one to construct a family of solutions that are connected with a particularly 
simple one by a finite chain of such transformations. 

The Backlund transformations play a central role in our method serving as an al- 
ternative to the standard Bethe ansatz technique. The nested Bethe ansatz solution of 
GL(iT)-invariant spin chains consists in successive increasing the rank of the group by 
applying the Bethe ansatz repeatedly. In this way, one can descend from GL(K) to 
GL{K — 1) until the problem gets trivialized at K = 0. At intermediate stages of this 
procedure, one introduces a number of auxiliary "Q-functions" (eigenvalues of Baxter's 
Q-operators) connected with the T-functions by Baxter's TQ-relations. Their zeros with 
respect to the spectral parameter obey the system of Bethe equations. This purely quan- 
tum technique has a remarkable classical interpretation [TJ [2] : it is equivalent to a chain 
of Backlund transformations for the Hirota equation. The TQ-relations appear then as a 
constituent of auxiliary linear problems for the Hirota equation. The rank of the group 
becomes an additional discrete variable, with the dependence on this variable being again 
described by the Hirota-like equation. Since the solutions are polynomials in the spectral 
parameter, their zeros obey equations of motion of a finite-dimensional dynamical system 
in discrete time. The equations of motion are just Bethe equations. 

Recently, this approach was applied [9] to supersymmetric spin chains constructed 
by means of GL(K |M)-invariant solutions to the graded Yang-Baxter equation [TOj, [TT] . 
In this case, there are two rather then one additional discrete flows corresponding to 
the bosonic and fermionic ranks, K and M. Their consistency leads to a non-trivial 



2 



bilinear relation for eigenvalues of Baxter's Q-operators (the QQ-relation). In the present 
paper we extend these results to GL(i^|M)-invariant spin chains with twisted (quasi- 
periodic) boundary conditions. The twisting parameters enter the solution as continuous 
parameters of the Backlund transformations. 

It is implied that the reader is familiar with the standard notions and facts related to 
supergroups and their representations [T2"j [T3| H%] , as well as with the quantum inverse 
scattering method constructions [151 EE fTTj . 

2 The TT-relation 

Let us recall the construction of the family of commuting transfer matrices in integrable 
spin chains. It is basically the same for ordinary and supersymmetric models. We con- 
sider lattice models (spin chains) with the symmetry supergroup GL(K\M) constructed 
by means of GL(i^|M)-invariant /^-matrices. Such R- matrices depend on the spectral 
parameter u G C and act in the tensor product Vq £g> V\ of two linear spaces, where 
irreducible representations tt and 7Ti of GL(K\M) are defined. The space Vq is usually 
called auxiliary space and V\ (local) quantum space. For our purposes we need the case 
when tci is the vector representation, i.e., V\ = V = C K © C M while ir is an arbitrary 
tensor representation of the supergroup. The i?-matrix reads 

R 01 (u) = u + 2j2(-l) m ME af s) ® e Pa . (1) 

a/3 

Here p(j3) is parity of the index f3 (p(/3) = or 1), e Qj g are matrices with the entries 
( e a/3)a'f3' = 8aa'8/3f3', and TtQ(E al3 ) are generators of the gl(K\M) superalgebra in the 
representation tcq. The first (scalar) term is to be understood as u multiplied by the 
unity matrix vro(J) ® where / G GL(K\M) is the unity element in the group and 
Iy x is the identity operator in the space V\. This i?-matrix is the GL(i^|M)-invariant 
solution to the graded Yang-Baxter equation [TO, [LT] • The supergroup invariance means 
that 

xo(g) ® m(g) Roi(u) = R i(u) n (g) ® tti(^) (2) 

for any g G GL(K\M). 

In order to introduce generalized integrable spin chains on N sites, take N copies of 
the space V = V\ = V 2 = . . . = Vn (one for each site of the chain) and the corresponding 
i?-matrices Raiiu) acting in V ® V{. The Hilbert space of states of the model, H, is 
the tensor product of local quantum spaces Vi over all sites of the chain: H = ®f = iVi. 
We will call 7i the quantum space of the model. The quantum monodromy matrix is 
constructed as the product of the i?-matrices Roi(u) in the space V : 

T(u) = R 01 (u - £i)R 02 (u - f 2 ) . . . R 0N (u - £ N ) . (3) 

The quantities £j are input data characterizing the (inhomogeneous) spin chain. The 
supertrace of the quantum monodromy matrix taken in the space Vq gives a family of 
operators in the quantum space depending on u and ir Q (called transfer matrices), which 
mutually commute for any values of these parameters. A more general construction 
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involves twisted (quasi-periodic) boundary conditions defined by means of a diagonal 
matrix 

g = diag {x 1 , . . . , x K , y ± , . . . , y M ) e GL(K\M) (4) 

(for simplicity, we identify elements of the supergroup GL(K\M) with matrices from its 
vector representation). Then the commuting family of transfer matrices is given by the 
formula 

T^\u;g) = stT m (Tr (g)T(u)) . (5) 

The graded Yang-Baxter equation combined with the GL(K\M) invariance implies that 
they commute for different u and tc (but not for different g\): [T^ n °\u; g), T^(u'\ g)\ = 
0. 

Below in this paper we will be especially interested in the case when the represen- 
tations 7To are "rectangular", i.e. correspond to rectangular Young diagram. Given a 
rectangular diagram of length s and height a, let 7r" be the corresponding representation. 
We define the quantum transfer matrices T(a, s, u) for rectangular representations in the 
auxiliary space by the formula 

T(a, s, u) = str^o (jc^(g)T(u — s+afj . (6) 

It differs from (jSJ) by the shift of u which is convenient in what follows. As a rule, we 
will not indicate the dependence on g explicitly. One may formally extend this definition 
to zero values of a and s which correspond to the trivial representation ttq = 7r° = 7T0 
(^(fiO = 1 f° r ari Y 9 ^ GL(K\M)). Taking into account that ^(E^) = for all 
generators of the superalgebra, we conclude from formulas (CQ), (jHJ) and (jSJ) that 

N N 

T(0, s, u) = n {u - s - Q , T(a, 0, u) = JJ (« + a - Q (7) 

i=l 3=1 

where ^ are the same as in (j3J). So we see that T(0, s, u) and T(a, 0, u) are unity operators 
in the quantum space multiplied by scalar polynomial functions. These functions are fixed 
input data of the problem. We put 

N 

0H = no-&)> ( 8 ) 

then T(0, s, u) = <p(u — s), T(a, 0, u) = <p(u + a). 

The transfer matrices constructed above are linearly independent but are connected 
by non-linear functional relations. In particular, the transfer matrices ([6]) for rectangular 
representations are known [51 [6] to obey the TT-relation 

T(a, s,u+l)T(a, s,u—l) = T(a, s+1, u)T(a, s — 1, u) + T(a+1, s, u)T(a—l, s, u). (9) 

This is the famous Hirota bilinear difference equation [8], where T plays the role of the r- 
function. Since all the transfer matrices mutually commute, they can be simultaneously 
diagonalized by a w-independent similarity transformation, and thus the same relation 
is valid for any of their eigenvalues. Keeping this in mind, we will think of the transfer 
matrices as scalar functions and call them T-functions. Our strategy is to treat equation 
dH]) as the basic equation of the quantum theory trying to derive the results for the 
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spectrum of quantum transfer matrices from it. To do this, we need to specify the 
boundary conditions and analytic properties of the solutions. The boundary conditions 
will be discussed in section 4. The analytic properties in u are determined by the type 
of the R- matrix (see a more detailed discussion in [3]). For quantum spin chains with 
polynomial i?-matrices (CQ) all the T-functions T(a,s,u) are polynomials in u. Thus we 
are led to the study of polynomial solutions to the Hirota equation. 

One should note that the general solution to the Hirota equation with the boundary T- 
functions at a = and s = as in ([7]) is given by the Bazhanov-Reshetikhin determinant 
formula [41. In our normalization it has the form 



where the function H(u,a) is defined as follows: H(u,0) = l/0(w), H(u,1) = 1, 
) = nj^i 1 <f)(u + a - 21) at a > 2. It can be expressed through the gamma-function 
for any a: 



Given arbitrary functions T(l, s, u), formula (Till gives a solution to the Hirota equation. 
However, in general this solution does not obey the required analytical and boundary 
conditions. In particular, the right hand side at a > 2 has apparent poles at zeros of 
the function H(u — s,a). Their cancelation is possible if the polynomials T(l,s,u) are 
chosen in a special way. 

We conclude this section by a few words on normalization of the solutions. It is easy 
to check that the transformation 

T(a,s,u) — > f (u + s + a)fi(u + s — a)f 2 (u — s + a)f 3 (u — s — a)T(a, s, u), (12) 

where /j are arbitrary functions, leaves the form of the Hirota equation unchanged. One 
may choose certain normalization of the solutions by fixing these functions in one or 
another way. In our normalization, all the polynomials T(a, s, u) have one and the same 
degree N equal to the number of sites in the spin chain. This formally includes special 
cases when one or more zeros of some of these polynomials are placed at infinity, then the 
degree is actually less then N. Other ways of the normalization are discussed in [3], [9]. 
In general position, the solutions (J5J) are irreducible, i.e., T(a, s, u) is not divisible by any 
polynomial of the form f (u+s+a)f 1 (u+s—a)f2(u—s+a)f 3 (u—s—a), where at least one 
of the polynomials fi(u) has degree greater than 0. 



3 The N = case: characters of the supergroup 



Before proceeding further, it is instructive to consider the case N = which appears to 
be rather meaningful, though simple, and thus provides useful analogies and motivations 
for dealing with more complicated models. It also emerges as the u — > oo limit of models 
with N > 0. 



T(a,s,u) = H (u—s,a) det T(l, s+i — j, u+a+l—i—j) 



(10) 




(11) 



GL(K\M) 
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At N = there are no spin degrees of freedom and the T-functions do not depend on u. 
As is seen from the definition, they coincide with characters of the element g G GL(K\M) 
for rectangular representations: T(a,s,u) = x{ a i s \g)- The characters depend on the 
parameters Xi,i/j entering the matrix g (j4]). We assume that all Xi,yj are distinct non- 
zero numbers. For brevity, we will write x(a, s) instead of x{ a ) s\g) when g is fixed. 

Let us introduce the rational function 



YC=i{i -y m t) 

n£=i(i-z fc t) 



w (0 = v5T 1 ; i (13) 



where t is an auxiliary variable. As it follows from the character formulas for supergroups 
[13], w(t) is the generating function of the characters x(l> s ) while the inverse function 
w; _1 (£) is the generating function of the characters x( a > 1) : 

oo oo 

w(t) = J2x(l,s)t s , w-\t) = Y / (-±Tx(aA)t a . (14) 

s=l a=l 

The other characters (super- analogs of Schur functions) are expressed through x(l, s) or 
x(a, 1) by the Jacobi-Trudi determinant formulas: 

x(°> s) = det s+i-j) = det x(a+i~3 , 1) • (15) 

They imply the bilinear relation for the characters of rectangular representations: 

X 2 (a, s) = x(a + 1, s)x(a -l,s) + x(a, s + l)x(a, s - 1) (16) 

which is the w-independent version of eq. Q). It is known as the discrete KdV equa- 
tion (see, e.g., [HJ EE]) written in the bilinear form. Let us also mention the integral 
representation of the characters 



l<j<k<a n=l 



where it is assumed that all singularities of the function w(t) are outside the unit circle. 
This a-fold integral coincides with the partition function of the asymmetric unitary matrix 
model written through the eigenvalues. 

One can formally extend the definition of characters to negative values of a, s by 
putting them equal to zero. Because x(0, s) = x( a i 0) = 1 at a, s > 0, this is consistent 
with equation (IIS]) everywhere except at the point a = s = 0. To make the equation 
valid in the whole (a, s) plane, one should put either x(0,n) = 1 or x(n, 0) = 1 for any 
n G Z, all other x( a > s ) with negative a or s being zero. We choose the former option 
(consistent with the integral representation (fTTj) ). Using representation (fl5l) . it is not 
difficult to show that if a > K + 1 and s > M + 1 then x(a, s) = 0. Summarizing, we 
have: 

X(a, s) = if : 

(18) 

(i) a < or (ii) a > and s < 0, or (iii) a > K and s > M. 

We see that the domain where x{ a i s ) do not vanish identically has the shape of a "fat 
hook" formed by the union of the half-strips a>0,0<s<M and s>0,0<a<K 
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together with the horizontal ray a = 0, s < from the origin to minus infinity (see 
Fig. [Q. We denote this domain by H(K,M). We call the boundaries at a = and 
s = 0, a > exterior and the boundaries inside the right upper quadrant interior ones. 

The characters on the interior boundaries can be found explicitly. A simple calculation 
shows that x{^i M + n) and x[K + n, M) at n > are given by the formulas 

/ k \ K M 

x(K,M + n) = [Ux n k )Y[Y[(x i -y j ), (19) 

\fc=l / i=l 7=1 



/ M \ K M 

X (k + n , m) = n (-i/m) n n n - 

\m=l / i=l j'=l 



(20) 



(in fact this is a particular case of a more general factorization property, see [20J, chapter 
I, section 3, example 23). Therefore, they are connected by the relation 

x(K,M + n) = (-ir M (sdetg) n x(K + n,M), n > 0, (21) 

where 

x x x 2 ...x K 

sdet g = . 

VW2 ■■■Vm 

Using the determinant identities (Pliicker relations) for minors of a rectangular ma- 
trix built from the elementary characters x(l> s )> one can prove the three-term bilinear 
relations 

X(a + 1, s)x(a, s) - x(a, s)x(a + 1, s) = zx(a + 1, s - l)x(a, s + 1) , 

(22) 

X(a, s + l)x(a, s) - x(a, s)x(a, s + 1) = zx{a + 1, s)x(a - 1, s + 1) 

between the characters x( a ; s ) with the generating function w(t) and the characters 
x(a, s) with the generating function w(t) = (1 — zt)w{t). Both x and x obey the discrete 
KdV equation (|T6|) and thus relations (122!) generate the Backlund transformation for it. 



7 



Their meaning is to relate characters of GL(/c|m)-representations with different k,m. 
Indeed, let g k ^ m e GL{k\m) be the diagonal matrix 

g Km = diag (x ± , . . . ,x k , y u . . . , y m ) (23) 

obtained from g = Qk,m (ED by removing K — k eigenvalues Xk, Xk-i, ■ ■ ■ , %k+i and M — m 
eigenvalues y M , Vm-i, • • • , Z/m+i- We put 

Xk,m(a, s) = x(a, s|flfc, m ) (24) 

and, in the same way as in ( fl3l) . (tHJ), introduce the generating function of the characters 
Xk,m(M): 

Wfcm (t)= n F a ~^ ■ (25) 

The obvious recurrence relations 

w k -i, m (t) = (1 - x fc t) w k>m (t) , 

(26) 

W fc , m+ i(t) = (1 - y m +l^) Wk,m(t) 

allow one to apply equations (|22|) . where one should put x( a ) s ) = Xfc,m(a, s), x( a > s ) = 
XJfe-i,m(a, s) at 2; = x fc or x(a, s) = Xk,m-i(a, s), x(a, s) = Xfc,m(«, s) at z = y m . 



4 The boundary conditions 

In general, the Hirota equation has many solutions of very different natures. The most 
important additional ingredient which selects the class of solutions relevant to quan- 
tum integrable models is the boundary conditions in the variables a, s. Qualitatively, 
these conditions are the same as those for characters of rectangular representations of 
supergroups discussed in the previous section, and can be derived by means of a similar 
reasoning. Furthermore, the explicit formula flTJ for the i?-matrix implies that the high- 
est coefficient of the polynomial T-function coincides with the corresponding character: 
T(a, s, u) = x( a , s ) uN + 0{u N ~ l ) as u — » oo. 

The distinctive feature of the solutions of our interest is that they are required to 
vanish identically in some parts of the (a, s)-plane. The domain where they do not do so 
depends on the symmetry algebra of the quantum model. For the GL(i^|M)-invariant 
super symmetric spin chains it is the domain H(K, M) introduced in the previous section 
(Fig.[T]). (For the GL(i^)-invariant spin chains it degenerates to the half-strip < a < K, 
s > together with the two rays s = 0, a > K and a = 0, s < 0.) Similarly to (1181) . we 
can write: 

T(a, s, u) — if : 

(27) 

(i) a < or (ii) a > and s < 0, or (iii) a > K and s > M. 

As is easy to see, the shape of H(i^, M) is consistent with the Hirota equation in the whole 
(a, s)-plane. Although only the points with non-negative a, s have the direct physical 
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interpretation, it is important to consider the full domain (the "fat hook" complemented 
by the ray) since otherwise the Hirota equation would brake down at the corner point at 
the origin. 

The boundary values of T-functions have a rather special factorized form fixed by 
consistency with the Hirota equation. Indeed, on the boundary one of the two terms in 
the right hand side of eq. (Q vanishes resulting in constraints on the boundary values. 
For example, at a = eq. (Q converts into 



which is a discrete version of the d'Alembert equation with the general solution T(0, s, u) = 
f+(u + s)/_(tt — s) where f± are arbitrary functions. In a similar way, T(a,0,u) = 
f+{u + a) f-{u — a). Equations © show that in our normalization f+(u) = f-(u) = 1 
while f-(u) = f+(u) = <f)(u), where <f>(u) is defined in (jHj). A similar factorization holds 
true on the interior boundaries. However, in general both functions in the product are 
non-trivial and are to be determined from the TT-relation. The case of the ordinary group 
GL(K) = GL(K\0), when the vertical interior boundary coincides with the exterior one, 
is special in this respect. In this case all the boundary values enter as fixed input data. 
In particular, on the horizontal interior boundary we have T(K, s, u) = (f)(u + s + K). 

On the interior boundaries, we impose the condition 



which is consistent with the corresponding condition (|2ip for characters of the supergroup. 
For the periodic case (when g is the unity matrix) this condition was pointed out in [TJ. 
Equality (|28|) means that the values of the T-functions at the points of the interior 
boundaries equally spaced from the corner point differ by a constant factor only. Note 
that this condition trivially holds in the case M = 0, where these values are fixed from 
the very beginning as input data. As it was already mentioned, at any K, M > they 
are no longer fixed but are to be determined together with the T-functions inside the 
domain H(K\M). 

Concluding this section, we note that the parameters Xk, y m do not enter explicitly 
neither the Hirota equation nor the boundary conditions for it. This means that the 
Hirota equation with the fixed boundary conditions on the exterior boundaries of the 
form ([7J has a continuous (K + M)-parametric family of polynomial solutions. 

5 Auxiliary linear problems 

Like almost all known nonlinear integrable equations, the Hirota equation is a compat- 
ibility condition for an over-determined system of linear problems [2TI [1] . To introduce 
them, it is convenient to pass to the "chiral" variables 



T(0, s,u+l)T(0,s,u-l) = T(0,s+l,u)T(0,s-l,u) 



T(K, M + n,u) = (-l) nM (sdet g) n T(K + n,M,u), n > 0, 



(28) 



p — \{u — s — a) 



q = \{u + s + a) 



(29) 



T = o 



\{—u - s + a) . 
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The original variables a, s, u will be referred to as "laboratory" ones. Here are the for- 
mulas for the inverse transformation, 

a = q + r , s = —p — r , u = p + q (30) 

and for the transformation of the vector fields: 

d p = d u -d s , d q = d u + d a , d r = d a -d s (31) 

We set r(p, q, r) = T(q + r, —p — r, p + q) and introduce the following linear problems 
for an auxiliary function ip = ip(p, q, r): 



V(r + 1) + z T ^)\ r ^ l) T ^ = ^p + 1) 
T{p + l)r(r + 1) 

^ + i)- ^ ( ^l;7 1 l> = ^ + 1 ) 

T{q + l)r(r + 1) 



(32) 



where z is a parameter. In these formulas, we indicate explicitly only those variables 
that are subject to shifts. Using these equations, the function -0(p + 1,9 + 1) can be 
represented as a linear combination of ip(r), ip(r + 1) and ip{r + 2) in two different 
ways. Compatibility of the linear problems means that the results are to be equal. 
Equating the two expressions we see that the terms proportional to z 2 ip(r) and ip(r + 2) 
cancel automatically while the terms proportional to zip(r + l) yield a non-trivial relation 
(provided ip(r + 1) does not vanish) 

r(p+l,r + 2) r(r + 1) r(p+l, q+1, r + l)r(p + 1) 
r(p+l,r+l)r(r + 2) r(p+l, q+l)(p+l, r+1) 

T(p+l,q+l,r+l)r{q + l) r(q+l,r+2) r(r + 1) 



T(p+l,q+l)(q+l,r+l) r(q+l, r+l)r(r + 2) ' 
This equality states that the function 

T(p + l)r(g + 1, r + 1) + r(g + l)r(p + 1, r + 1) 
r(r + l)r(p + l,g + 1) 

is a periodic function of r with period 1 and an arbitrary function of p, q. Because no 
special periodicity is implied, we set this function to be r-independent. Therefore, we 
arrive at the relation 

r(p + l)r(q + 1, r + 1) + r(q + l)r(p + 1, r + 1) = h(2p, 2g)r(r + l)r(p + 1, q + 1) , 

where /i can be an arbitrary function of p and g. In the original variables this equation 
reads 

T(a + l)T(a - 1) +T(s+ l)T(s - 1) = h{u-s-a, u + s + a) T{u + l)T(u - 1) . 

From the boundary conditions ([7]) at a = or s = it follows that h = 1 and we 
obtain the Hirota equation ([9]). Note that the parameter z entering the linear problems 
disappears from the non-linear equation. In fact this is clear from the very beginning 
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because z can be eliminated from equations fl32l) by the transformation ip — > z p+q+r if). 
Nevertheless, we keep this parameter because it will be important in what follows. 

An advantage of the "crural" variables is their separation in the linear problems: the 
first problem does not involve q while the second one does not involve p. However, in 
contrast to the "laboratory" variables a, s, u, they have no immediate physical meaning. 
Coming back to the "laboratory" variables, we set ip(p,q,r) = + r, — p — r, p + q) 
and rewrite the linear problems fl32|) in the form 



(33) 



t/ x T(a — 1, s+l,u)T(a, s — l,u + l) T . , , N T , ^ lN 

(a, s, u) + z^—^ ' 7 * a- 1, s+l, u = #(a-l, s, w+1 

i {a,s,u) 1 (a— 1, s, «+lj 

t/ x T(a — 1, s+1, u)T(a+l, s, w+1) _ . , , . _ , „. 

* a, s, u) - z V ' ' / \ ' ; * a- 1, s+l, u) = * a, s+l, m+I . 

i [a,s,u) I [a, s+l, w+lj 

In general, compatibility of linear problems implies existence of a continuous family of 
common solutions. As we have seen, the structure of the coefficient functions in our 
case is such that the compatibility is equivalent to the existence of at least one common 
solution (see [22], where this fact was pointed out in another context). 

Because the T-functions can vanish identically at some a, s, we eliminate the denom- 
inators by passing to the new auxiliary function F = T^>, in terms of which we have 

T(a— l, s, u+l)F(a, s, u) + zT(a, s — l, u+l)F(a—l, s+l, u) = T(a, s, u)F(a— l, s, w+1) 

T(a, s+l, w+l)F(a, s, w) — zT(a+l, s, u+l)F(a—l, s+l, u) = T(a, s, u)F(a, s+l, u+1) . 

(34) 

Note that the second equation can be obtained from the first one by the transforma- 
tion T(a,s,u) — ► (— l) as T(— s, — a, u) (and the same for F) which leaves the Hirota 
equation invariant. However, the Hirota equation written for the function T(a,s,u) = 
(— 1) 2 < y°- 2+s2 )T{a, s,u) is form-invariant with respect to a larger symmetry group consist- 
ing of any permutations and changing signs of the variables a, s, u but the system of the 
linear problems (1341) is not. In fact the symmetry is realized in an implicit way. To make 
it explicit, we write the pair of equations (|34"1) in a matrix form, 



T(a—l,s,u) zT(a,s — l,u) \ / F(a, s, u — 1) 
T(a,s+l,u) —zT(a+l,s,u) J \ F(a— 1, s+l, u— 1) 

/ F(a — l,s,u) 



(35) 



T(a, s,u—l] 



V F(a,s+l,u) 



and multiply both sides by the matrix inverse to the one in the left hand side. Using the 
TT-relation, we get another pair of linear problems, 

T(a+1, s+l, u)F(a, s, u) — zT(a+l, s, u+l)F(a, s+l, u—1) = T(a, s, u)F(a+l, s+l, u) 
T(a, s, u+l)F(a, s, u—1) — T(a, s—l, u) F(a, s+l, u) = T(a+1, s, u)F(a—l, s, u) 

(36) 
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T(a, s, u) 



T(a, s, w) 



-T(a, s + 1, tt) 



T(a — 1, s, u) 







T(a+l,s,w) ' 
T(a, s — 1, w) 
— T(a, s, m + 

/ 



(37) 



which are equivalent to (and thus compatible with) the pair (|34p by construction. The set 
of four linear problems ([31]), fl3"ET) possesses the required symmetry. The Hirota equation 
can be derived as a compatibility condition for any two linear problems of these four, 
and the other two hold automatically. The four linear equations can be combined into a 
single matrix equation 

( F(a-l,s,u) \ 

F(a, s+1, u) 

F(a, s, it— 1) 

y zF(a—l, s+1, u-1) J 
where T(o, s, u) is the antisymmetric matrix 
/ T(a,s,u— 1) 

—T(a,s,u—1) 

T(a,s+l,u) —T(a—l,s,u 

i — T(a+l,s,u) —T(a,s — l,u) T(a,s,u+1 



(3f 



The Hirota equation implies that its determinant vanishes and rank of this matrix equals 
2. The symmetric form of the linear problems for the Hirota equation was suggested in 
j. For more details on the linear problems and their symmetries see [2"T| [T51 12"3"t |2"1] . 



6 Backlund transformations 

There is a remarkable duality between T(a, s, u) and F(a, s, u) [2TJ [1]: one can exchange 
their roles and treat eqs. (|34p as an over-determined system of linear problems for the 
function T with coefficients F. Their compatibility condition is the same Hirota equation 
for F: 

F(a, s,u+l)F(a, s,u-l) = F(a, s + l,u)F(a, s-l,u) + F(a + l, s,u)F(a-l, s,u) . (39) 

We thus conclude that any solution to the linear problems (I3"4"|) . where the T-function 
obeys the Hirota equation, provides an (auto) Backlund transformation, i.e., a transfor- 
mation that sends a solution of the nonlinear integrable equation to another solution of 
the same equation. 

Let us rewrite the linear problems (134"1) changing the order of the terms and shifting 
the variables: 

T(a+1, s, u)F(a, s, u+1) — T(a, s, u+l)F(a+l, s, u) = zT(a+l, s—l,u+i)F(a, s+1, u) 

T(a, s+1, u+l)F(a, s, u) — T(a, s, u)F(a, s+1, u+1) = zT(a+l, s, u+l)F(a— 1, s+1, u). 

(40) 
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T(u) 



F(u) 



T(u+1L 



F(u+1) 



T(u+1) 



T(u) 



T(u+1L 



F(u) 



^T(u+1) 



T(u) 



Figure 2: The graphical representation of equations HW in the (a, s) -lattice. Here a and 
s coordinates correspond to the vertical and horizontal axis respectively. 



These equations are graphically represented in Fig. [2] in the (a, s)-plane. They constitute 
the Backlund transformation T — > F in the bilinear form [HI [2TJ- Given a family of 
polynomials T(a, s, u) obeying the Hirota equation, one may pose the problem of finding 
polynomial solutions to equations (H0|) . It is easy to see that these equations are not 
compatible with the boundary conditions for F(a,s,u) and T(a,s,u) of the "fat hook" 
type with the same K and M. Indeed, applying these equations in the corner point of 
the interior boundary, one sees that if K, M for T and F are the same, then the boundary 
values must vanish identically However, it is straightforward to verify that equations 
10]) are compatible with the boundary conditions of the following two types: 



F(a,s,u) = if: 

(i) a < or (ii) a > and s < 0, or (iii) a > K — l and s > M, 



(41) 



or 



F(a,s,u) = if 



(i) a < or (ii) a > and s < 0, or (iii) a > K and s > M + 1 . 



(42) 



They are again of the "fat hook" type but with the shifts K—tK—lorM—tM + l. We 
refer to the corresponding transformations as BT^~ and BT^": F(a, s, u) = BT^~ (T(a, s, u)) 
for (T4T]) and F(a, s, u) = BT^~(T(a, s, u)) for fT421 . They depend on the parameter z. The 
same formulas fT40l) define inverse transformations BT^ = (BT]") _1 and BT^ = (BT^) _1 
if one treats them as linear equations for T with given F . In a more explicit way, the 
transformations BT^, BT^ are defined by formulas fj471) . fj48|) below. 

Repeating these transformations several times, we arrive at the hierarchy of functions 
Tk, m (a, s,u) (k = 0, 1, . . . , K, m = 0, 1, . . . , M) such that: 



a) They obey the Hirota equation in a, s, u for any k, m; 
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b) They are connected by the Backhand transformations 

T k -i,m{a,s,u) = BT^ (T k>m (a,s,u)), 

(43) 

T k ,m-i{a, s, u) = BT 2 (T k>m (a, s,u)); 

c) At k = K, m = M we have 

T KM (a, s, u) = T(a, s, u) . 

For the irreducible polynomial solutions this list should be supplemented by the "initial 
condition" To t o(a, s,u) = 1 at a = or at s = 0, a > and T 0t0 (a, s,u) = otherwise. 
The levels of the hierarchy are labeled by the pair of numbers k, m. The lowest level is 
0, while the highest one is K, M. We shall say that the T-functions T kiTn (a, s, u) belong 
to the level k, m. 

It is important to note that the parameter z involved in the definition of the Backlund 
transformations can be different for transformations of the two types introduced above 
as well as for successive transformations of the same type. We choose these parameters 
to be eigenvalues of the matrix g: z = x k for the transition (k,m) — > (k — l,m) and 
z = y m for (k,m) — > (k, m — 1). More precisely, the transformation BT^ is written in 
the form 

T k ,m{a + 1, s, u)T k - 1>m (a, s,u + l) - T ktTn (a, s,u+ l)T k _ hm (a + 1, s, u) 
= %k T kt7n {a + l,s-l,u+ l)T k -. 1>m (a, s + l,u), 

(44) 

T kjm (a, s + l,u + l)T fc _i )m (a, s, u) - T kjTn (a, s, u)T k _ hm (a, s + 1, u + 1) 
= x k T ktTn (a + l,s,u+ l)T k _ ljm (a - 1, s + 1, u), 
while BTg in the form 

T fcj7n _i(a + 1, s, u)T k>m (a, s, u + 1) - T fcjm _i(a, s, u + l)Tfc, m (o + 1, s, w) 

= y m T k , m -i{a + 1, s - 1, u + l)T fcjm (a, s + l,u), 

(45) 

7fc,m-i(a, s + 1, u + l)T k>m (a, s, u) - T^^a, s, u)T k>m (a, s + 1, u + 1) 

= y m Tk, m -\{a + l,s,u+ l)T k>m (a — l,s + l,u) 

(k = 1, . . . ,K, m = 1, . . . , M). If one ignores the w-dependence (which disappears in 
the u — > oo limit), then these formulas become the bilinear relations between characters 
mentioned in Section 3 (see equations ( l22|) and comments after them). It is easy to notice 
that each of the equations in (144p . ( 1431) is actually a dynamical equation for a function 
of three variables rather than five. For example, the first equation in ( 1441) acts in the 
subspaces m = const and u + s + a = const. Upon restriction to the corresponding 
three-dimensional hyperplanes in the linear space with coordinates a,s,u,k,m, each of 
these equations can be put in the standard Hirota form by a linear change of variables. 
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For completeness, let us give a symmetric description of the Backlund transformations 
through the matrix equations of the form (j57|) . For the direct and inverse transformations 
we need to introduce two antisymmetric 4x4 matrices T < - ±1 - ) (a, s, u) of the type 



T (e) (a,s,w) 



/ T(a, s, u — e) —T(a,s+e,u) T(a+e,s,u) \ 

—T(a,s,u — e) T(a — e, s,u) T(a,s — e,u) 

T(a,s+e,u) —T(a — e,s,u) — T(a, s,u+e 

i — T(a+e, s,«) — T(a, s — e, it) T(a, s,u+e) y 



(46) 

where e = ±1. Then the transformations T k _i tTn = BT 1 (T& m ), T k ^ m+ i = BT^"(Tfc m ) are 
defined by the matrix equation 



/ Tk-i, m {a-l, s, u) 



k,m+ 



i(a— 1, s, w) ^ 



(47) 



7fc_i, m (o, s+1, it) T fc)m+ i(a, s+1, u) 

T fc _i jm (a, s, u-1) T fe)m+ i(a, s, u— 1) 

^ x fc T fc _ ljm (a-l, u-1) y m+1 T k}jn+1 (a-l, s+1, u-l) J 

and the inverse transformations T fe+l m = BTf(T k>m ), T k>m _i = BTg (T fc m ) are defined 
by the equation 



/ 7fe + i, m (a+l, s, w) 
T fc+lj7n (a, s-1, w) 
7fc+i, m (a, s, u+1) 



7fc,m-l(a + l, s, w) \ 

Tfc, m _i(o, s-1, w) 
T fcim _i(a, s, w+1) 



(4? 



y x fc+ iT fe+ljm (a+l, s-1, w+1) j/ m T fc)m _i(oH-l,s-l,tH-l) / 

where T^(a, s,u) are the matrices f|46l) with entries at the level k, m. 

Moreover, a careful analysis of the equations (14*4")) - (148]) shows that boundary values 



of the T-functions at each level k, m are subject to the same relations as at the highest 
level. On the exterior boundaries, the T-functions have the specific form similar to (J7J), 
i.e., T k )m (0, s, u) is a function of u — s while T fc>m (a, 0, u) is the same function of u + a. 
Let us introduce the special notation for them: 



T fcj7n (0, s,u) = Q k ,m{u 



Tu m (a,0,u) 



Q 



k,m 



U + a) . 



(49) 



The polynomials Qk, m {u) play a very important role. They will be identified with eigen- 
values of Baxter's Q-operators. The polynomial Qk,m( u ) = 4>( u ) ls a fixed input data 
which determines the model. The polynomials Qk,m{u) at lower levels are to be found 
in the course of solution. In analogy with (jSJ), we fix their highest coefficients to be 1. 
Applying (I44p . (1451) on the interior boundaries, we conclude that if (128)) is valid at the 
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highest level, then T-functions on the interior boundaries at each level k, m are connected 
by a similar relation: 

T k , m (k,m + n,u) = (-l) nm (sdet g k , m ) n T k , m (k + n } m } u) } n > 0. (50) 

The matrix g kjTn e GL{k\m) is defined in (j2"5|) and sdet g k>m — X\ . . . x k j (yi - ■ ■ y m ). 

We see now that by applying BT^~ and BT2 to a solution of the Hirota equation in 
the domain H(.K"|M) with the boundary conditions ( 1271) . we can successively transform 
it to the trivial solution in the degenerate domain H(0|0). This "undressing procedure" 
allows one to construct solutions to the original problem, as it will be shown below. 



7 Recurrence relations for the operator generating 
series 

The Backlund transformations from the previous section can be reformulated in the 
operator form as recurrence relations for difference operators of infinite order. Let us 
consider the following operator: 

m^g /' 1 ^;/- 1 ^ (51) 

(the common denominator is introduced for normalization). It serves as a non-commutative 
generating series for the T-functions T(l, s, u). Similar objects can be introduced at any 
level k,m: 



m, m (u) = E Tfc ' m(1 n g, V' 1} e 2s9 " • (52) 

It is the generating series for the functions T fc m (l, s, u). Clearly, Wo,o( M ) = 1 (recall that 
T 0:0 (l, s, u) — unless s = and T (l, 0, u) = Qo,o(u + 1) = 1). We also note that the 
series formally inverse to fl5"2l generates the T-functions T fcm (a, l,u): 

m» = E(-ir e 2aa - Tk ' m{ ^ X Z ~_ a ~ 1} ■ (53) 

a>0 tyk,m\ u Z J 

For the proof, see [9]. 
Set 

v , , Qk, m (u+2)Q 

X Km [u) = x k — —— — — , (54) 

,r , S Qk,mr-l{U + 2) Q km (U-2) 

Y k> m{u) = y m — ~ r —^ l — f-s — • (55) 

Qk,m-l{U) Qk,m[ U ) 

Using the linear problems at a = 0, it is a straightforward calculation to prove the 
following recurrence relations for the operators W kym {u): 

W k -i, m (u) = (l - X Km (u)e 2d -) W k , m (u) , 

(56) 

W k , m+ i(u) = (l - Y Km+l {u)e 2d -) W k , m (u) . 
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These formulas are operator (and u-dependent) analogs of fl26|) . The shift operator 
t = e 29u plays the role of the variable t while the variable u is absent in ( 1261) . Notice also 
that Xk, m (u), Yk tm (u) turn into Xk, y m in the limit u — > oo. 

Here we present some details of the proof (see also [9J, where a slightly different 
version of the recurrence relations is proved). Consider the first relation. We have: 



w fc -i, m (u) - yv k ,m(u) = Y, 



s>0 



7fe-l,m(l, S, U + S - 1) Tfc m (l, S, M + S - 1) 



Qk-l,m(u) 



Qk,m(u) 



J2.sd u 



• (57) 



To transform the expression in the square brackets, we rewrite the first equation in (JIH 
at a = in the form 

7fc,m(l, S, It + S - 1) Tfc_i jm (l, s,u + s-l) 



Xk 



Qk,m(u) Qk-l,rn(u) 

Q k ,m{u + 2)Q k 2) T k , m (l,s-l,u+s) 



Qk,m(u)Qk-l,m(u) 



Qk,m(u + 2) 



and continue the equality: 



Because Tj^ m (l, — !,«) = 0, the sum in right hand side can be written in the form 



(5f 



T kjJn (l,s-l,u+s) 2sdu _ 2du ^ T fe , m (l,s,M + s- 1) 



s>0 Qk,m\U' 2) s >q Qk,m{^ 



and the first recurrence relation is proved. The proof of the second one is completely 
similar. 



8 Factorization formulas and TQ-relations 

The recurrence relations established in the previous section allow one to represent the 
operator generating series (15 ip in a closed factorized form, where each factor contains the 
Q-functions only. Namely, Wk,m{u) can be obtained as a result of successive application 
of the recurrence relations fl56l) to Wo,o(u) = 1. In this way, moving first in the re- 
direction from (0,0) to (0, M) and then in the fc-direction from (0, M) to (K,M), we 
get: 

wk,m(u)= n (i-XkM^e^y 1 ■ n (i-vw^) (59) 

K>k>l M>m>l 

where the ordered product is defined as 

fj A = Aj Aj_! . . . A I+1 Aj . 

J>i>I 
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Applying the recurrence relations in the different order (fc-direction first, m-direction 
next), we arrive at a different but equivalent representation: 



II (l-^,m(w)e 29 »)- II (l-X kfi (u)e 29 -y . (60) 




In fact one can apply the relations fl56l) in any other order determined by a chosen 
zigzag path from the point (0,0) to the point (K,M). Each step in the fc-direction, 



(k,m) — > {k + l,m), brings the factor (1 — X k+ i tm {u)e 29u J while each step in the re- 



direction, (k,m) — > (k,m + 1), brings the factor ( 1 — Y k , m+ i{u)e 2du \ which are to be 



multiplied according to the order of the steps. This yields many other ways to factorize 
the operator Wk,m (m). Their equivalence follows from the compatibility of the recurrence 
relations fl56l) which is expressed by the discrete "zero curvature" condition 



(l - Y k ^ m+1 (u)e 29 -) (l - X k , m (u)e 29 ") = (l - X k , n+1 (u)e 29 ") (l - Y Km+l (u)e 2d -) 



on the (k, m)-lattice. The two sides of this equality correspond to two different ways to 
obtain W k -i, m +x{u) from W k , m {u). 

The equalities fl59l) and fl60l) as well as the similar equalities with different orderings 
are generalized Baxter's TQ-relations in a generating form. Equating coefficients in 
front of different powers of the operator e 2du , one obtains expressions for the T-functions 
T(l,s,w) = Tx,Af(l, s, u) through X kjm , Y k>m and thus through the Q-functions Q k , m (u) 
with 1 < k < K, 1 < m < M. For example, the simplest TQ-relation contained in (!59j) 
has the form 



( |55|) . The zeros of the latter are to be constrained by the system of Bethe equations 
derived in the next section. 

9 QQ-relation and Bethe equations 

Our starting point in this section is the discrete zero curvature condition (|6T|) . Comparing 
coefficients in front of different powers of the shift operator, we note that 

Yk~i,m+i{u)X ktm (u + 2) = X ktm+1 (u)Y kjm+1 (u + 2) 

holds identically, and thus get the only non-trivial relation 



which after the substitution (15^1) . fl55l) becomes a functional equation for the Q- functions. 
As a simple calculation shows, it is equivalent to the following bilinear equation (the 





(61) 




Yk-l,m+l( u ) + X k>m (u) — X k ^ m+ i(u) + Y k ^ m+ i(u) , 



QQ-relation" [9]): 



XkQk-l,m-l(u)Qk,m(u + 2 ) - y m Qk,m(u)Qk-l,m-l(u + 2) 



(63) 



) Qk-l,m(u)Qk,m-l(u + 2) 
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for the polynomial functions 

QtAu) = II (« - u f' m) ) ■ (64) 
i=i 

Let us remark that the QQ-relation (|63|) can be recast to the standard form of the 
Hirota bilinear difference equation in "crural" variables u, m, —k by passing to the func- 
tion 

Qk,m(u) = a fc , m e (/3fc+ ^> Q k , m (u) , (65) 
where the new parameters f3k and 7 m are related to the Xk, y m by the formulas 

Xk = e Wk-fo-l) j y m = e 2(7m-l-Tta) _ ( 66 ) 

They are fixed uniquely by putting /3 = 7o = 0. This transformation eliminates the 
coefficients x k , y m in (as well as in (154]) . (1531) ) and, with a proper choice of Ofc )m , the 
QQ-relation acquires the coefficient-free form 

Qk-l,m-l(u)Qk,m(u + 2) ~ Qk,m(u)Qk-l,ra-l(u + 2) = Qfc-i, m (tt)Qfc, m _i(tl+ 2) (67) 

suggested in [9]. 

The QQ-relation fl63|) provides the easiest and the most transparent way to derive 
Bethe equations for roots of the polynomials Qk,m(u). Putting u in (1631) successively 
equal to the roots of each Q-function entering the equation, one obtains a number of 
relations which, after some rearranging, can be written in the form 

gfc-ig {uf m) ) Q Km [uf m) - 2) Q k+ltm (uf m) + 2) 
Q k -i,m \uf m) - 2) Q Km {uf m) + 2) Q k+1>m (uf m) ) 



Xk+1 



(6* 



Q k ,m+l (uf'^) Qk,m (Uj' - 2) Qk,m-1 (uf + 2 
Qk,m+1 (uf ,m) - 2) Q Km (uf ,m) + 2) Qfc, m _i (uf' m) 



Vm+1 



(69) 



^ / (k,m)\ ^ ( (k,m) \ 

Vfc+l,m I J l°2k,m-l {Uj + ZJ Xk+1 

Qk+l,m (uf' m) + 2) Q k , m -i {uf^) Vm 



(70) 



Qk,m+1 (uf^) Qk-l,m (uf'^ ~ 2) = y m+1 

Qk, m +i (uf' m) - 2) Q fc _i ;m (uf ' m) ) ^fc 

They hold inside the K x M rectangle in the (k, m)-lattice and serve as elementary 
building blocks for systems of Bethe equations. Each such system corresponds to a 
zigzag "undressing" path from (K, M) to (0, 0). On the parts of the path (k + 1, m) — ► 
(A;, m) — > (A; — 1, m), (k, m + 1) — > (A;, m) — > (A;, m — 1), 1, m,) — * (A;, m) — ► (k,m — 1) 
and (A;, m + 1) — ► (k,m) — > (A; — 1, m) one uses (1681) . (1691) . flTOj) and ([71]) respectively. 
These systems are different but equivalent. For a more detailed discussion on this point, 
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see [9]. As an example, we give here the chain of the Bethe equations for the simplest path 
(K, M) — > (0, M) — > (0, 0). Moving from (K, M) to (0, M), we have the equations 

(72) 



Qk-iM (4' M) ) Qk,M {uf M) - 2) Q k+ltM (uf^ + 2 

Q k -i, M \uf M) - 2) Q KM (uf M) + 2) Q k+lM (uf M) ) ' ' x k+1 



where k — 1, . . . , K — 1. They agree with the chain of Bethe equations presented in 
for the bosonic case. At the turning point, the equation is 

Vm ' 

Finally, moving from (0, M) to (0, 0), we have the equations 

Qo,m+l ( u f' m) ) Qo,m [Uj - 2) Q ,m-1 (uj + 2) y m+1 



Qi, M (4°' M) ) 


! Qo,M-l 1 


(«r> + 2) 


Qi,m ( u f' M) 


+ 2) 


! <5o,A/-l 1 





(73) 



Qo,m+l (V-°' m) - 2) Q ,m (w5°' m) + 2) Qo,m-l («j 

where m = 1, 2, . . . , M — 1. 



(0,m) 
i 



(74) 



10 Conclusion 

We have obtained a solution of the TT-relation (the Hirota equation) obeying all the 
required boundary and analytic conditions. The solution is given by the determinant 
formula (1101) . where the polynomials T(l,s, u) are determined by the expansion (1511) 
of the factorized operator (|59p . The coefficients of the latter are expressed through 
the Q- functions via equations fl54l) . fl55|) . where the roots of the polynomials Qk,m(u) 
are constrained by Bethe equations which result from the bilinear QQ- relation fl63l) . It 
should be noted that the solution is not unique. Given boundary conditions, there is a 
finite set of solutions corresponding to different quantum states of the generalized spin 
chain. 

We emphasize that this solution, typical for quantum problems solvable by Bethe 
ansatz, has been obtained by purely classical methods of the theory of soliton equations on 
the lattice. The key role is played by Backhand transformations for the Hirota difference 
equation. From the viewpoint of the theory of classical soliton equations, our method 
consists in constructing a chain of successive Backlund transformations which reduces 
the problem to a trivial one. Each transformation of this chain involves a continuous 
parameter (the classical spectral parameter) which is identified with an eigenvalue of the 
matrix which defines the twisted boundary conditions in the quantum integrable model. 
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